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Abstract 

We present a prescription to calculate manifestly gauge invariant tree-level he- 
licity amplitudes for arbitrary scattering processes with off-shell initial-state gluons 
within the kinematics of high-energy scattering. We show that it is equivalent to Li- 
patov's effective action approach, and show its computational potential through nu- 



merical calculations for scattering processes with several particles in the final state. 



1 Introduction 



Besides determining the precise nature of the recently discovered Higgs-like boson, the scien- 
tific plans of the Large Hadron Collider (LHC) include further tests of complex dynamics of 
the Standard Model and searches of physics beyond the Standard Model d). Such a scientific 
program requires also theoretical control over background and signal processes occuring at the 
collisions. A substantial part of the necessary calculations involves the application of Quantum 
Chromodynamics (QCD). It enables calculations for LHC physics via various factorization the- 
orems, allowing for the decomposition of a given process into a part characterizing the colliding 
hadrons, called the parton densities, and a so-called hard part characterizing the parton parton 
scattering and production process. The predictive power of QCD largely relies on the fact that 
this hard part can be calculated within perturbation theory. Regarding the factorization theo- 
rems, we will focus on high-energy factorization [2] here, which applies when the energy scale 
involved in the scattering process is high and larger then any other scale in the process. For recent 
applications see [f3]-[9]] . This approach is particularly interesting as far as LHC physics is con- 
cerned since the energy is indeed the largest available scale and therefore this approach should 
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provide the optimal method to calculate cross sections for various processes. Furthermore, high- 
energy factorization already at lowest perturbatively defined order provides kinematical effects 
in hard matrix elements which are relegated in other approaches to corrections of higher orders. 
The evolution equations of high energy factorization providing transversal momentum dependent 
parton densities sum up logarithms of energy accompanied by a coupling constant lfTO]4T3l . Be- 
cause of the dependence on the parton's transversal momentum, parton densities of high-energy 
factorization have to be convoluted with the hard process which is calculated with the initiating 
gluons being off-shell. 

Several frameworks exist to calculate gauge invariant matrix elements for partonic processes 
initiated by off-shell gluons lfl4l4P71 . In order for such a framework to be widely used, a practical 
numerical tool would be desirable, like the ones that already exist for the calculation of tree- 
level on-shell scattering amplitudes, e.g. |fT8l - [22l . They operate on the amplitude level and use 
helicity methods, and thus are very efficient and universal, and can deal with essentially an 
arbitrary number of external gluons. The present study is a generalization of ll23ll . and provides 
the proof of concept for the application of the helicity method to high-energy factorization with 
two initial state gluons being off-shell. It includes an implementation of this method in a form 
of a numerical tool. In particular we develop a new prescription for calculating gauge invariant 
high-energy factorizable amplitudes for any tree-level process with two gluons in the initial state 
being off-shell. Our results turn out to be in agreement with results obtained using Lipatov's 
action. 

In short, our prescription is established as follows. It is known []24| that the usual Feynman 
graphs contributing to an amplitude in an on-shell calculation are in general not sufficient to 
obtain gauge invariant results if any of the external legs are off-shell. Thus the off-shell pro- 
cess needs to be embedded into a larger on-shell process, from which a gauge invariant off-shell 
amplitude needs to be disentangled. The first crucial point of our approach is that we manage 
to have exactly the kinematics of high-energy factorization for the off-shell process, and at the 
same time have exact on-shellness for the embedding, through the application of complex mo- 
menta ll25l for the on-shell "external" partons of the embedding. The second crucial point is that 
we have a well-defined prescription how to decouple the unneeded on-shell partons. 

The paper is organized as follows. In the Section |2] we introduce the kinematics of high- 
energy factorization and introduce our prescription which provides gauge invariant amplitudes 
for any tree-level process with off-shell initial-state gluons. In Section [3] we provide the deriva- 
tion of our prescription using the spinor helicity formalism and analytic continuation of the am- 
plitude to a region where the momenta of the auxiliary partons are complex. Furthermore, we 
show that the obtained amplitudes obey the correct collinear limit when the initial- state gluons 
become on-shell. In Section |4] we discuss the connection of our method with our former re- 
sult [|23l for g*g — > Ng. In Section [5] we show that our prescription reproduces the results 
obtained from the effective action approach, i.e. that the effective vertices of [fT51 correspond to 
our amplitudes with reggeized gluons. In Section |6] we present the numerical application of our 
framework to various LHC processes with a simple model for the transversal-momentum depen- 
dent parton density function. We conclude in Section [7] and the appendices contain details of the 
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spinor formalism and explicit calculations of reggeon-gluon-gluon and reggeon-reggeon-gluon 
vertices. 



2 Prescription 

We start with stating the main result of this paper, which is a prescription to calculate scattering 
amplitudes in quasi-multi-Regge kinematics in the high-energy limit, i.e. scattering amplitudes 
with off-shell initial-state gluons, and with no kinematical restrictions on the final-state particles 
other than momentum conservation and on-shellness. The derivation will be given in the next 
section. 

The exact definition of the scattering amplitudes we deal with requires the introduction of 
two momenta t-\ , l 2 such that 

l v ii=l 2 .l 2 = , £ r £ 2 =2E 2 , "E>0. (1) 

The high-energy limit is enforced by demanding that the momenta of the off-shell initial-state 
gluons are given by 

k^ = x^ + k^ , ]% = x 2 t$ + k$ ± (2) 



for xi,X2 G [0, 1] and 



«i -lci± = «2-ku = «1 -TciL = ^2-k 2 ± = . (3) 



The expressions for the off-shell momenta are restricted in the sense that for a component 
proportional to £ 2 is missing, and for k 2 a component proportional to is missing. This is the 
high-energy limit, or high-energy approximation. In our derivation, however, it will not be the 
result of an approximation. The kinematics of Eq. © is set from the start. The prescription to 
calculate the scattering amplitude for the process 

g*(ki)g*(Tc 2 ) -> X(k, + k 2 = p, +p 2 + . . . + p n ) 

with n particles in the final state is to embed the process into a quark-scattering process, and 
apply eikonal Feynman rules to the quark lines. More precisely: 



1. Consider the process qA q B — > qA qB X, where qA, qB are dis- 
tinguishable massless quarks not occurring in X, and with mo- 
mentum flow as if the momenta Pa, Pb of the initial-state quarks 
and Pa', Pb' of the final-state quarks are given by 

Pa — K i ' Pb — K 2 > Pa /— Pb , — u - 



_ _Pa' 
1 



2 
n 

Pb^ " Pb' 




Associate the number 1 instead of spinors with the end points of the A-quark line, 
interpret every vertex on the A-quark line as g s T^ instead of — ig s Ty 
interpret every propagator on the A-quark line as 5y/£i -p instead of iS^/p'. 
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3. Associate the number 1 instead of spinors with the end points of the B-quark line, 
interpret every vertex on the B-quark line as g s Ty instead of — ig s T^ y H , 
interpret every propagator on the B-quark line as 5tj instead of ibi)/i>. 

4. Multiply the amplitude with 

p = ixiy/-2k^ ^ ix 2V /-2k| ± ^ 

For the rest, normal Feynman rules apply. The coupling constant in the denominators of © are 
obviously necessary because the quark- scattering process carries a higher power of the coupling 
constant. The factor includes y/l x y/l to correct for the difference in color representation 
between the quarks and gluons, and a factor ixi yf— kf x ix 2 yj — k| * s necessary to obtain the 
correct collinear limit. 

3 Derivation 

The prescription involves embedding the process with the off-shell initial-state gluons into the 
quark- scattering process qAqB — > 9a9bX, with eikonal Feynman rules for the quark-lines, and 
deformed kinematics. The starting point of the derivation will be the scattering amplitude for the 
on-shell process, with normal quark lines 

A{c\ K c^ -» qAqBX) . (5) 

This is a well-defined gauge invariant function of the external momenta, helicities and colors. 
The key-point of our derivation is that we will not apply any approximations througout, and 
will manifestly keep gauge invariance throughout. The high-energy limit of Eq. © will be 
set from the start. We do need, however, to extrapolate the amplitude © beyond its physical 
interpretation. More explicitly, we will use the fact that in order to guarantee gauge invariance 
of © the external momenta need to be on-shell and need to satisfy momentum conservation, but 
are not required to be real (cf. I1251 ). 

3.1 Gauge invariant helicity amplitudes 

Rather than putting everything in an appendix, we state some essential points regarding the for- 
malism here. Consider the following decomposition of four- momenta [|26l . Given two light-like 
momenta l\ , li such that l\-lz ^ 0, we define 

i\ = \UY\U\ , l\ = \{U\y*U. (6) 

The exact definition of the spinors in the expressions above are given in Appendix [A] We have 
U-U= h-h = iris = U-U = U -h =lrU = h-h = h-U = , (7) 
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Figure 1: A few terms in the classification of the graphs contributing to q A qB — > c|a X w.r.t. 
the gluons attached to the quark lines. 



while 

l 3 -U = -lvk- (8) 
Any four- vector p can now be decomposed as 

^-^■P f , , ivy fV . U-y ^ try p (9 , 

V ~ i r i 2 1 «r«2 2 Mi 3 Mi 4 ' 

The first two terms are the usual terms in the Sudakov parametrization in terms of t-\ , li, and the 
last two terms form an explicit decomposition of the transversal component. If t-\ , li are real, 
then the momenta £3, l\ are in general complex. The spinors in Appendix [A] are well-defined also 
for complex momenta, and we have the following identities 



<*3l 


= (ti\ 


, l«3l 


= W 




= t«ii 


, l«3> 


= l<2> 


{U\ 


= («ii 


, U 


= |e 2 ] 


[U\ 


= [h\ 


, l«4> 


= Ki> 



(10) 



The amplitude of the processs q A qB — > qA qB X contains Feynman graphs that consist of a 
quark-line of type A with a single gluon attached and a quark-line of type B with a single gluon 
attached. They are represented by the first term on the r.h.s. in Fig. [U The heuristic picture is 
that these gluons are going to play the role of the off-shell gluons in g*g* — > X. It needs to be 
stressed that we will not drop any of the other graphs and that we keep the full gauge invariant 
sum. We demand that 

Pa-Pa'=K =xi^ + k^ , pg-pS,=lcjf = x 2 «J + lc£ L . (11) 

We do not set p A equal to £1 . We are interested in obtaining a scattering amplitude for the 
process g*g* — > X, and are not concerned about the momenta of the quarks. We only need to 
make sure we have overall momentum conservation and on-shellness. At first sight, the latter 
requirement seems to be in contradiction with requirements (ITTb . and it would be if we required 
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the momenta of the quarks to be real. Choosing, however, 

p£ = (A + x 2 )^+K 4 2 ^ , v »,=All-4q (12) 



with 



K 1 — on ' K l 



M 2 ' «1'«2 



3 _ _M<2± 4 _ _ V]<2± n . 

K2_ ' K2_ ' (13) 



both on-shellness of Pa>Pa'>Pb>Pb' and Eq. (fTTj) are assured to hold, for any value of A. Fur- 
thermore, following the identities (flOT ). we have 



<p A 'l oc , |p A ] oc 

(Pb'|oc(£ 2 | , |p B ]oc|£ 2 ], (14) 

so we may simply assign the following spinors to the quarks in the amplitude without spoiling 
gauge invariance: 

qA(PA) -> |£il , qA(pA') -M^il 

q B (p B )^|£ 2 ] , q B (pB')^(*2l. (15) 

As a result, in the graphs of the first term on the r.h.s. of Fig. [H both quark-lines are connected 
to gluons with eikonal vertices 

(ii\y»\l,]=2i» , (£ 2 |yHy=2^. (16) 

Other graphs contain A-quark propagators and/or B-quark propagators, and depend on A. We 
want to stress that the constructed amplitude is gauge invariant for any value of A. 
Now, we consider the limit 

A — > 00 , (17) 

and we will see that it is a well-defined limit. Considering the role of A in Eq. (fT2l) . one might 
interprete it as the high-energy limit. Realize, however, that A is a dimensionless parameter, not 
the energy, and that the limit is not an approximation. The actual high-energy approximation 
was essentially made by requiring Eq. © from the start. The limit (TT7T ) should rather be seen as 
the final step to extract a physical amplitude, independent of imaginary momentum components. 
The momentum p of an A-quark propagator can be written as 

p^ = (A + x p )^ + V 2 +P±> (18) 

for some numbers x p ,t) p and (complex) pj_, which are independent of A. Since l-\ -p A ' = 0, we 
can write 

y v l v l 2 = £1 -p = U ■ (p - PaO = «i V , (19) 
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where p' = p — p A ' is the momentum that would be flowing in the propagator if the external 
quarks had the momenta p A = k] and p A / = 0. The square of the momentum is given by 

p 2 = 2(A + x p )-y p M 2 + pi, (20) 
so for the propagator of an A-quark line with momentum p we find 

A— >co 

p 2 " " 2(A + Xp)y p « r « 2 + pi 2y p £ r «2 2£ r p' " 

Similarly, the propagator of a B-quark line with momentum p will become / 2 /(2^2 - p') 5 where 
p' = p — p B ' is the momentum that would be flowing in the propagator if the external quarks 
had the momenta Pb = k 2 and ps' = 0. Any gluon attached to an A-or B-quark line will do so 
via an eikonal vertex, e.g. for the A-quark line we have 

(I, | y * J** u- ■ ■ %\ = («i I y * («i I Y 1 UkU I • • • %\ 

= (2C)(2^)-.- . (22) 

The factors (— 2i) from the vertices and (i/2) from the propagators on the eikonal lines cancel, 
except for one vertex on each line. This happens in each graph, and thus leads to an overall factor 
2x2, which we remove. This concludes the derivation of the prescription up to the last point 
about the extra factor. 



3.2 Matching to the collinear limit 

Let us denote the obtained amplitude by 



A 



T-A)A> l B)B 



iqAqB -> qAqBX) 



(23) 



The sub-scripts indicate the color degrees of freedom of the quarks. As explained earlier, the 
coupling constant in the denominators is obviously necessary in Eq. ©. We will now show that 
the amplitude also has to be multiplied with the numerator. In order to arrive at a satisfactory 
helicity amplitude for the process g*g* — > X, the obtained amplitude needs to be matched at the 
collinear limit of kjj_, k| ± — > to the process g g — > X with on-shell initial-state gluons. We 
perform the matching at the level of the squared amplitude, averaged over initial- state helicities 
for the on-shell process, and averaged over the azimuthal angle of the transversal components 
of the off-shell gluon momenta for the off-shell process [|27l . The color content of the obtained 
amplitude is different from the content of the on-shell process, and the matching of these requires 
also summation over initial-state colors. So the factor F in Eq. © should be such that 



1 



a u a 2 Ai,A2 



l-A,)A> l B,)B 



lim TA 



VUaAbJb 



;qAqB -> qAqBX) 



'1 ,<t>2 



(24) 
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Figure 2: An example of an off-shell current. The enumerated external lines are on-shell, the 
other gluon is off- shell with momentum k = pi + P2 + P3- 

The average over the azimuthal angles can be defined e.g. for k 1± by parameterizing it in terms 
of |k] j_| and 4^ following 

k ^ = ^|=^( ei * l£ 3+ e " i *^4) , (25) 
and performing the integral over cbi , dividing the result by 2n. 

3.2.1 Color 

The off-shell gluons in the process g*g* — > X carry color degrees of freedom in the adjoint 
representation. The additional quarks in the process c|a c| b - > q A q B Xarein the fundamental 
representation. Eventually, the amplitude needs to be squared and summed over color degrees 
of freedom. Abreviating A a ' h {g*g* — > X) just to A a ' h , where a, b are the color indices of the 
off-shell gluons, we have 

Y_ \A a *\ 2 = y_ A ^ 5QC5bd ( AC4 Y 

a,b a,b,c,d 

= Y_ Aa,h (2TrT Q T c ) (2TrT b T d ) (yi c - d )* = Y_ | 2 A'j,ki| 2 , (26) 

a,b,c,d i,j,lc,l 

where 

A' i]M (g*g* X) = ^T^T k b v ^ b (g*g* X) . (27) 

a,b 

So we see that attaching the quark lines to the off-shell gluons does not alter the color content 
of the amplitude apart from a factor 2 if, eventually, we sum over the color degrees of freedom 
of the quarks. As an alternative approach, we could simply interpret the amplitude to be in the 
color- flow representation [12811291 . which would also require an extra factor yl x yl. This 
explains the factor y/2 x \fl in Eq. (|4]). 

3.2.2 Collinear limit 

We consider the amplitude of Eq. (|23l . Let JJJ be an off- shell current, i.e. the sum of all connected 
sub-graphs that contain a given set of the external particles plus one off-shell gluon. An example 



8 



is depicted in Fig. [2] The propagator of the off-shell gluon is not included in the off-shell current. 
The momentum k H of this gluon is the sum of the momenta of the external particles in J" Gauge 
invariance implies current conservation, which dictates that 

k^ = 0. (28) 

If this off-shell current is attached to the A-quark line via the numerator of a gluon propagator 
d^(k, n) in an axial gauge with gauge vector rO 1 , we have 

—(TO + ^KJS). (29) 

For the gluon with momentum k^ 1 = x\l^ + k^ ± , we have l\ -ki = 0, and we see that the corre- 
sponding off-shell current is always attached to the A-quark line as (— ^J^) , independently 
from the gauge. Using current conservation again, k^J^ — 0, we have 

- «i% = ^ K Jl • (30) 

Including the remaining piece of propagator i/k 2 ± , we see that the off- shell gluon with momen- 
tum ki is attached to the A-quark line via the contraction with 

— i k^ 

(3D 



The second part in the expression above will give the correct average over the azimuthal angle, 
and the first part represents the reciprocal contribution to the factor of Eq. ©. Similarly, for the 
gluon with momentum k^ = X2^ + we ^ n< ^ mat ^ * s attached to the B-quark line via the 
contraction with 

— i k^ 

- Kl± (32) 



W- k L V-H± 

For kj ± , k| ± — > 0, contributions that do not contain both off- shell gluons with the momenta 
ki , k2 vanish, and only those from the first term on the r.h.s. of Fig. [Qremain, which are exactly 
the ones one would take into account for the on-shell process gg — > X. The only difference is 
that these gluons are contracted with (13TT ) and (|32l instead of polarization vectors. We consider 
the collinear limit at the level of the squared amplitude, averaged over the azimuthal angles of 
ki_i_ and k 2 ^. Notice that for k^ ± , k 2± — > 0, the amplitude only depends on ku_, k 2 j_ through the 
contractions with (l3TT ) and (l32l) . and it is enough to prove that the contractions of the squared 
amplitude with 

'*I.iAl\ /^2±\ (33) 



k 2 / ' \ -k 2 



V.1 
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lead to the correct average over initial-state polarizations. We just consider ku_, the other one 
goes analogously. Using its decomposition < fT3l in terms of £3, £4, we have 



V* 1 V y i^ 3 P^P y k- 4 P^P y P^P y -4- P^P y 

-k 2 ± ~ K?2£r£ 2 K ?2£ r £ 2 2£ r £ 2 * 

Rembering Eq. ©, we see that £3/^1 '^2 and -it are just the two polarization vectors 

in the Kleiss-Stirling [30] construction for light-like momentum £1 with reference momentum 
£2, and that the third term on the r.h.s. in Eq. (|34l l gives the desired polarization average. The 
dependence of k,/ic| on 4) is proportional to e 21 *, so that the first two terms vanish in the average 
over (b. 



4 Processes with one off-shell gluon 

For processes with only one off- shell gluon, i.e. g* y — > X, with y = g, q, q, one has to consider 
the process q A ij — > q A X, and only apply the eikonal Feynman rules to the A-quark line. In this 
case, the helicity amplitude can be simplified further by using a special axial gauge. 

In the previous section we saw already that the off-shell current ]^ with momentum ki = 
Xi£i + ki^ is attached to the A-quark line as (— £^J^), independently from the gauge. Eq. (l29l) 
also tells us that if we choose the axial gauge with n 11 = l^, then all contributions to the amplitude 
with any other off-shell current attached to the A-quark line vanish. In other words, the only 
graphs with gluons attached to the A-quark line that need to be taken into account are one of the 
following: 

1. the only gluon attached to the A-quark line is the gluon with momentum k 1 , 1 = xi l\ + k^ ± 

2. all gluons attached to the A-quark line are on- shell 

Consequently, if X in the process g* y — > X does not consist of gluons only, the second option 
above cannot occur, and we see that 

• for the process g* y — > X where X does not consist of gluons only, helicity amplitudes 
can be calculated by applying Feynman rules without attaching a propagator, in the axial 
gauge with gauge vector £] , and by contracting the off-shell gluon with 

k^ 

gLL . (35) 

We assume the off-shell gluon to be in the adjoint color representation here. If X consists of 
gluons only, then the contribution coming from graphs with all on-shell gluons attached to the A- 
quark line is exactly the 'gauge-restoring' amplitude of equation (25) in [|23l . It can be eliminated 
by choosing the reference momentum of the polarization vector £^ of at least one of the on-shell 
gluons equal to l^, so that for this on-shell gluon i-\ ■ e = 0. 
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ki +p 2 




k, +p 2 P ^995> 



Pa' 



Figure 3: Feynman graphs for q A g — > q A g. 



A 



-A' 
n 



perm. 



O 
o 
o 
o 

J 

n 



A' 



Figure 4: Effective vertices consisting of n gluons attached to an eikonal line. The sum is over 
all permutations of (1 , 2, . . . , n). 



5 Effective reggeon vertices 

In this section we will show that the prescription for the calculation of scattering amplitudes with 
off-shell gluons can reproduce the results of [fl"51 . The effective vertices involving reggeized 
gluons presented in that paper correspond to amplitudes with off-shell gluons as presented here, 
in the case that the non-reggeized gluons are on-shell. Reggeized gluons in these vertices corre- 
spond to gluons with momenta restricted to be of the form of Eq. ©. To be more specific, in the 
following, we show that our prescription provides 

^&>fa§ x V2T£ x F^^ + Oc i; pi,p 2 ,...,Pn;fe) , 06) 

V _k l V _|C 2 

where V is the reggeon-reggeon-n-particle vertex from [fl5l . Analogously, the prescription pro- 
vides the single-reggeon vertices. 

We present two explicit calculations in Appendix El namely for g*g — > g and g*g* — > 
g. The necessary graphs in the embedding of g*g — > g are depicted in Fig. |3] The last two 
graphs in Fig. |3] contribute as an effective vertex of the type of Fig. HI These vertices contain all 
possible ways a number of gluons can be attached directly to the eikonal line. According to our 
prescription, p A / = 0, and the denominator factor of each propagator is given by — i-\ -p where p 
is the sum of the incoming momenta of the gluons to the right of the propagator. We have 

A - „ 'A' 




where we abreviate 

£i = «rPi, (38) 
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and where 



jai JQ2 . . . JQn 

G a ,a 2 -a n (^l,^2,-..,^n) = M /_ , f , , ; , Z Z r , (39) 

- c,llc,l + c,2j • ■ • + t,2 + r C, n -1 J 



perm. 



which is identical to equation (19) in [[T5l . The sum is over all permutations of (1 , 2, . . . , n). 
Due to momentum conservation and i-\ -ki = 0, the arguments to G QlQ2 ... Qn are restricted such 
that £,! + £,2 H + £, n -i + £, n = 0. In AppendixOwe prove that this implies the identity 

Tr{G aiai ... a j£, 1 ,£, 2 ,...,£ >n )} = 0, (40) 

so that 

(Ga,.^^,^,...,^ =2T^Tr{T c G QlQl ... an (^,£, 2 ,...,£, n )} . (41) 
Combining the foregoing, we see that 



km? n 




" (-igr 1 ) A:£g::z(iut2,...,u) , (42) 



with 

2k 2 

A-X-a^(^i^2, • • ^n) = — ^ W • • ■ C Tr{T c G Qiai ... an (^i, £, 2) • • • , U) , (43) 

which is the equivalent of equation (18) in |fT5ll for the eikonal (— )-direction. The factor E in the 
denominator is necessary to compensate for the fact that the vector l-\ has the dimension of E, 
while the vector nr in lfT5l is dimensionless. For n = 3, the factor ig 2 in Eq. (|42|) corresponds 
to the same factor in the denominator on the l.h.s. of equation (43) in |[T5l . The remaining minus 
sign is due to the fact that in the formula above, all momenta are considered to be incoming, 
while pi , p 2 in equation (43) in |fT5l are outgoing. Equation (41) in |fT5l contains a misprint, and 
the factor — i should be accompanying the As. 

The vertices of Fig. @] contain all eikonal vertices and propagators. Any vertex involving 
reggeized gluons can be constructed using these in combination with normal Feynman rules. 
By recursively following the graphical decompositions of Fig. [5] to Fig. [H it can be understood 
that our prescription delivers the reggeon-reggeon-particles vertices up to 3 particles. On one 
hand, by expanding all blobs and effective vertices, one can see that all Feynman graphs for the 
process qA qB — > qA qB X are obtained. On the other hand, one can easily identify the graphs in 
our figures with those in figures 4, 5 and 8 from [fl~5l . Only for our Fig. [8] our notation is slightly 
more compact, and the comparison requires the expansion of the first term on the second line 
following Fig.fTJ 
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Figure 5: Graphs for the reggeon-gluon-gluon vertex. 




Figure 6: Graphs for the reggeon-reggeon-gluon vertex. 
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Figure 7: Graphs for the reggeon-reggeon-gluon-gluon vertex. 



13 




Figure 8: Graphs for the reggeon-reggeon-gluon-gluon-gluon vertex. 



6 Numerical application 



We implemented the prescription of Section |2] into a numerical program to calculate differential 
cross sections within high-energy factorization for arbitrary final states: 



do 



dxi 



d 2 k 1± f(x b k 1± ] 



dx 7 



d 2 k 2± f (x 2 , k 2± ) d<D (ki , k 2 ;pi , p 2 , . . . , - 



>n) 



(N 



1 



L L 

colors helicities 



\A(g*(k,)g*(k z )^X[ VhV2 ,...,pn))\ 
2sxfx 2 Sym(X) 



(44) 



where k] , k 2 are defined as in Eq. © and the sums are over all final-state helicities, and all 
colors, including the initial-state colors. Sym(X) is the symmetry factor associated with possible 
identical particles in the final state, and s is the square of the total center-of-mass energy. The 
function f is the transversal-momentum dependent gluon density function. 

Tree-level helicity amplitudes A are automatically calculated following the recursive Dyson- 
Schwinger method, in the spirit of 111 8 11201121 113 Ul32l . The helicity and color sums are performed 
numerically, evaluating so-called color-dressed helicity amplitudes Il3~3~ll3"4~l . As mentioned at 
the beginning, this happens following the prescription of Section [2l i.e., the recursive Dyson- 
Schwinger relations are augmented with the eikonal Feynman rules, and amplitudes with two 
initial-state and two final-state eikonal quarks are evaluated. This increases the multiplicity by 
2, but this causes no serious burden for tree-level processes up to 6 final-state particles. Also, 
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the eikonal vertices and propagators are computationally less involved than the other vertices 
and propagators. The phase space integration of the final-state momenta is performed with 
KALEU I1351 , and the integration over the initial-state variables Xi , X2, k.u_, k 2 ^ is performed with 
adaptive grids following the method of 11361 . 

Since the scope of this paper is rather a proof of concept than a phenomenological study, we 
choose to present results within a simple model for the transversal-momentum dependent pdf. 
We use 



where Qo is a scale of a few GeV, and f co i (x, \i) is the collinear gluon density from CTEQ6L1 [I371 . 
For simplicity, we just fix the scale u. to the mass of the Z-boson. Eq. (1451) has the property that 
the off-shellness of the initial-state gluon is squeezed to zero for small values of Q , such that 



This way, we can check the collinear limit without having to worry about the collinear limit of 
the transversal-momentum dependent pdf. 

In the following, we present some differential distributions for the processes 




(45) 



Qo^O 



lim d 2 k ± f (x, kjj = x f co i(x, ]i) . 



(46) 



g*g* — > bbZ — > bb u+u. 
g*g* ->bbZg -> bb>V~ 9 
g*g* ^bbg 
g*g* -> bbgg . 



(47) 
(48) 
(49) 
(50) 



The Standard Model input parameters used are 



m w = 80.41 9 GeV r w = 2.12GeV 
m z = 91 .1 882 GeV T z = 2.4952 GeV 
= 1.16639 x 10 5 GeV 2 
a s = 0.13 



(51) 
(52) 
(53) 
(54) 



The value of the fine- structure constant is defined following the scheme 



sin 2 9w = 1 — Tnw/ m z 
V2G 



(55) 



(56) 



71 



The masses of the fermions involved in the processes are 



m b = 4.62 GeV 
m u = 0.1056583 GeV 



(57) 
(58) 
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All calculations are performed at a center-of-mass energy of 14TeV. The phase space of gluons 
in the final state is restricted such that p T > 1 GeV, and such that 

x/Acb 2 + Ay 1 > OA , (59) 

where Acj) is the difference in azimuthal angle and Ay the difference in rapidity with any of the 
other final- state gluons or quarks. 

In Fig. [9l Fig. [TOl Fig. [TT] and Fig. [12] some distributions are depicted for the mentioned 
processes. Results for two values of the scale Qo in Eq. (l45l) are presented: Qo = 0, which refers 
to the result within collinear factorization, and Qo = 4 GeV. The latter value was chosen such 
that it leads to a ratio of the cross sections for the process g*g* — > bbZ similar to those presented 
in 11381 . The absolute values are considerably smaller, because we present results for an off-shell 
Z-boson with only the p+pr-channel taken into account. 

We want to stress that the presented results should merely be seen as a proof of the compu- 
tational potential of our approach. Still, we may state that the results are consistent with those 
presented in e.g. 113811391 . The cross sections within high-energy factorization are suppressed 
compared to collinear factorization at leading-order QCD, and the transversal-momentum distri- 
butions show a larger spread. 

7 Summary 

We presented a prescription to calculate tree-level helicity amplitudes for arbitrary scattering 
processes with off-shell initial-state gluons within the kinematics of high-energy scattering. We 
showed that the amplitudes are manifestly gauge invariant, and that the approach is equivalent to 
Lipatov's effective action approach. We established the connection with earlier work regarding 
multi-gluon amplitudes with one off-shell initial-state gluon, and we showed the computational 
potential of the presented approach with numerical calculations for scattering processes with 
several particles in the final state. 
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p T (bb)[GeV] y(bb) 

Figure 9: Differential cross sections of the transverse momentum p T and the rapidity y of the Z- 
boson and the bb-pair in the process g*g* — > bbZ — > bb|j. + |jr. The different curves correspond 
to the different values of Qo in Eq. (1451) . with Q = refering to collinear factorization. The 
total cross sections are 26.181 ± 0.014pb for Q = OGeV, 16.159 ± 0.036 pb for Q = 4GeV. 
These numbers are considerably smaller than those in [38], because here the Z-boson is off-shell 
with only the |o. + pr -channel taken into account. 
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Figure 10: Differential cross sections of the transverse momentum p T and the rapidity y of the 
Z-boson and the bb-pair in the process g*g* — > bbZg — > bb(x + |j.-g. The different curves cor- 
respond to the different values of Qo in Eq. (1431) . with Qo = refering to collinear factorization. 
The total cross sections are 93.57 ± 0.21 pb for Q = GeV, 76.98 ± 0.23 pb for Q = 4 GeV. 
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Figure 11: Differential cross sections of the transverse momentum p T and the rapidity y of the 
bb-pair in the process g*g* — > bbg. The different curves correspond to the different values 
of Qo in Eq. (1451) . with Q = refering to collinear factorization. The total cross sections are 
695.41 ± 0.32 pb for Q = GeV, 1 71 .74 ± 0.20 pb for Q = 4 GeV. 
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Figure 12: Differential cross sections of the transverse momentum pj and the rapidity y of the 
bb-pair in the process g*g* — > bbgg. The different curves correspond to the different values 
of Qo in Eq. (|45l) . with Q = refering to collinear factorization. The total cross sections are 
772.4 ± 1 .2 |xb for Q = GeV, 393.9 ± 1 .2 jxb for Q = 4 GeV. 
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A Definition of spinors 

We introduce the bi- spinor of a four- momentum p following 

/ p - P3 -P1 + 1P2^ /aY . 

cr-p = o po - cr-p = r \ , > (60) 

\— Pi - 1P2 P0 + P3 / 

where o"o is the identity 2 x 2-matrix and a = (<7i , 02, 03) are the Pauli matrices. We also define 

&-p = o po + cr-p • (61) 
For p 2 = Pq — pf — p| — pf = 0, we may write 

o--p = L(p)Up) T , &-p = R(p)R*(p)\ (62) 
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with 



up) = ; ( " Pl + iVl ) uv) - v/|po+P3 ' f ~ V] " ip2% 



VlPo +Psl V Vo+Vs J Vo + Vs \ P0 + P3 / 

vSn(Po+Ps\ R, (p ) = ^^(P° + M . (63 ) 

Po + Ps VPi + m) VWfpsI \Pi - W 

If p is real with p > 0, then L*(p) = L(p)* and R*(p) = R(p)*, and L(p), R(p) are the well- 
known Weyl spinors. The spinors above are defined also for complex momentum, as long as it is 
on-shell. Introducing the bra-ket-notation, we write 



IP] 
IP) 



( U P) ) , <Pl = (o,U P ) T ) 

(r(p)) ' [PI = ( R *^ T >°) ( 64 ) 



lp)[pl + lp](pl=(.° 7)=^. (65) 



so that 

IrAPnl 4- MM = I 

&- v 

This leads to the Weyl representation of the y-algebra, with 

75 = IyVyV = (~* J) • (66) 

We state some basic relations. We trivially have 

[p|p) = (p|p] = , (67) 
and explicit calculation reveals that also 

(pip) = [p|p] = . (68) 
From explicit calculation also follows that 

(p|q)[q|p]=2p-q. (69) 

One more important relation is that 

2p.q = (p|q)[q|p] + <p|q](q|p] 
= <p|(|q)[q| + lq](ql)lp] 

= (v\i\v]- (70) 

By decomposing any arbitrary four-vector into two light-like vector, it can easily be shown that 
this relation also holds if q is not light-like. 

Realize, finally, that in explicit numerical calculations, it is usually not wise to identify pi = 
Px>P2 = Py>P3 = Pz> since in collider physics one may typically encounter situations in which 
po = — p z - Therefore, it is better to identify e.g. pi = p z , P2 = p x > P3 = Vy ■ 
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B Explict calculation of reggeon-gluon vertices 



In this appendix, we calculate some simple reggeon-gluon vertices using our prescription. The 
calculations are performed in the Feynman gauge. 



B.l The reggeon-gluon-gluon vertex 

First, consider the process g*g — > g. The process is embedded in 

<U(PA,j) 9(P2, a,v) -> q A (p A ',i) g(ki + p 2 ,b,a) . 

The arguments a and b refer to the color of the gluons, the arguments v and cr are the Lorentz 
indices of the gluons, and the arguments j and i refer to the color of the quarks. The necessary 
Feynman graphs are depicted in Fig. [3] According to the prescription, the fermion line should 
be interpreted as an eikonal line with momenta Pa = ki and p A ' = 0. The sum C of all three 
graphs is 



C = ig s x, a /-2k 2 



T**f fcab V ^ (kbp2) _ ki _ p2] 
K 1 



+ 



where V^ ya is the usual vertex for three vector-bosons 

v^( Pl , P2 , p 3 ) = (pi - V iY ^ y + iv2-vsV g w + (ps - pi V g^ • (71) 

Using £1 - ki = 0, we find 



C = ig s X! 



c .ccab 



[2(«, -p 2 )g w - (2k! + p 2 )X + (k, - P2 )T] 



«rp 2 



£|£? 



Inserting [T a , T b ] = if cab T c , we get 
-g s xi V2TI f cab 



C 



2(£i -p 2 )g w - (2k! + Vl yq + (k, - Vl yi 



k 2 
1 £|£? 



«rP2 



2|= y™" (P2, ^ k! + p 2 , b; k, , c) , 



-k 2 



with 



y]\ (p 2 ,a;pi,b;k b c) 



2(£, -p 2 ) g w - (2 Pl - Vl yi° - (2 P2 - Pl yq 



1c? 



«rp 2 



£|£? 



(72) 



which, for £i = E n , is identical to equation (34) in 02 
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B.2 The reggeon-reggeon-gluon vertex 

Next, we consider the process g*g* — > g. It is embedded in 

qA(PA,j) qB(pB,l) -> qA(pA',i) qB(pB',k) g(ki +k 2 ,b,o) , 

where the arguments j, I, i and k refer to the color of the quarks, the argument b refers to the 
color of the gluon, and a is the Lorentz index of the gluon. The necessary graphs are depicted in 
Fig. Ell According to our prescription, we have p A = ki, p B = k 2 , j>a> = Pb' =0, and the sum 
of the graphs is equal to 



C = -2g s xix 2 vkfkf 



(T a T b ) tj £T£f -i JQ , (T b T a )ij -i TQ 



£ l -(k 1 +k 2 ) Y\ 



£ r (- k2 ) k 2 T K Q l £ ^ 

i (T c T b ) kl £ 2 ^ , „, -i (T b T c ) k i« 



+ T S V k? Mk 1+ k 2 



kf «2-(-ici: 



Using £-| -ki = £ 2 -k 2 = 0, we find 



C = -2g s X!X 2 Wkfk 



TCTa xcab 
" Mj 'kl T 



[2(«! -k 2 )£f - 2(£ 2 -k, )£? + (£, -£ 2 ) (k, - k 2 r 



£vk. 



k2 T kl ^2v + Tjj i-\ ^ ^ 2 



i [T c ,T b ] kl ^ 



£ 2 -k! 



Inserting [T a , T b ] = if abc T c and [T c , T b ] = if cbQ T Q , we get 



C = -2g s xix 2 wkfk| 



TCTa -pcab 
" 'tj 'kl T 



[2(i r k 2 )q - 2{i 2 -h)^ + [i, -h) (k, - k 2 r] 



(£ 2 -k0k 
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Using the anti-symmetry of f abc , we finally find 



_ 2g s xix 2 l r l 2 bQ 

c — - - - Hj 'kl T 



VW 2 
XiX 2 E 2 



"Mi 



k 2 1 
W 1 £ 2 -k! 2 



k 2 



2T^T k a i r^(k 1 ,c;k 1 +k 2 ,b;k 2 ,a), 



with 



2g s f 



abc 



E 2 ' €rk 2 



«rk 2 



E 2 £ 2 -k! 



(k 2 



kil CT 



(73) 



which, for £] = E n and £ 2 = E n + , is identical to the first line of equation (41) in [fT51 combined 
with equation (42) in 03). Realize that £1 -£ 2 = 2E 2 . 



C Proof of Eq. <gQD 



Here, we proof that for a cycle-symmetric tensor T tl " ln and the set of the complex numbers 
£,1 , . . . , £, n satisfying £,1 + . . . + £, n = the following identity holds 

Y 7 r = , (74) 

{Wn} ^ ^1 + ^) • • • + h 2 + ■ • • + O 

where the sum goes over all the permutations of the indices i] , . . . , i n . 
Let us decompose the sum over all permutations as follows 

q-M ...in 

y r- = 2l Ti, '" ln X. Q ^1' •• -'T-n) , 

ft lv ..,w ^ ^ + ^2) • • • + ^ + • • • + £,i n _, ) , n} , _ . n) 

(75) 

where ^ 

a (ii, . . ..in) = -. r- (76) 

£, il (£, ii +^ i J...(£, i]+ £, i2 + ... + £, in _ ] ) 

and , . . . , i n }' denotes the set of all non-cyclic permutations of the indices. The notation 
(ii , . . . , in) refers to the set of all cyclic permutations. Let us construct the following mero- 
morphic function 

f (£) = Z ^— / r • (77) 

(-£,) [-£, + l h ) (-£, + ^ + h 2 ) . . . (-£, + ^ + h 2 + ■ . ■ + £,i n _, ) 
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The single proper poles are 

z = 0, 

Zn = + h 2 + ■ • ■ + h n _, ■ (78) 

The residues of f (£,) in those points are 

res f (£,) = — — — — 1 — r- = a(ii,...,i a ), (79) 

i=z a £,i, (£,i, + U 2 J . . . (£,i, +U 2 + ... + £, in _ 1 j 



res f f £.) = - 

(-**, ) ^ ft* + W ■ ■ • + • • • + ^ ) 

1 



h 2 (h 2 + ti 3 )... {h 2 + .-. + ) + • • • + h n _, + 

= a (i 2 , . . . , in, ii ) , (80) 



res f (£,) 



^2 {-h, ~ h 2 ) l~h 2 ) h 3 lh 3 + *m) • ■ ■ {h 2 + ... + O 

1 

= a (t3, . . . , i n , il , X2 ) , 

and so on. Let us now consider the integral of f (£,) over the circle with origin £, = and radius 
R — > oo. Obviously 

o f(£,)d£, = 0. (81) 



On the other hand we have 



= ^f(£,)= Y_ a(i b ...,i n ) (82) 



and thus the sums in Eq. (1751) vanish. Let us finally remark that we cross-checked the identity 
explicitly up to n = 6 with the help of FORM BOll . 
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